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Abstract 

The neutrino mass hierarchy can be determined, in principle, by measuring a phase in the 
disappearance oscillation probability in vacuum, without relying on the matter effect, using a single 
channel. This phase is not the same for the normal and inverted neutrino mass spectra. In this 
paper, we give a complete description and physics understanding of the method. The key feature 
of the method is to detect advancement (normal) or retardation (inverted) of the phase of the 
atmospheric-scale oscillation relative to the solar-scale oscillation. We then show that this method 
can be realized with the recently proposed resonant v e absorption reaction enhanced by Mossbauer 
effect. The unique feature of this setup is the ultra-monochromaticity of the observed v e 's. Because 
of this feature the phase advancement or retardation of atmospheric-scale neutrino oscillation is 
detectable after 20 or more oscillations if the source and the target are made sufficiently compact 
in size. A quantitative estimate of the sensitivity to mass hierarchy resolution is given. We have 
also examined how a possible continuation of such experiment can be carried out in order also to 
achieve high precision (few %) determination of the solar-scale oscillation parameters Am^ and 
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I. INTRODUCTION 



Among the remaining unknowns in neutrino physics, resolving the neutrino mass hierar- 
chy may be one of the key importance to our understanding of physics of neutrino masses 
and the lepton flavor mixing Given the significance of the problem, it is natural that 
many methods are proposed for determining the mass hierarchy. The most conventional one 
would be to do neutrino and antineutrino comparison in long-baseline (LBL) accelerator 
based oscillation experiments 0, [3[ . Because of the difference in how the matter effect inter- 
feres with the vacuum oscillation effect it can signal the sign of Am^, and hence the mass 
hierarchy. Other less conventional ideas include the methods using atmospheric neutrinos 
0] , v e and Up disappearance channels @, 0] , supernova neutrinos [§1, ^ , and neutrino- less 
double beta decay [l0(] . These lists of references cannot be complete, and some of the other 
relevant references are quoted therein. 

In this paper, we discuss yet another way of determining the neutrino mass hierarchy by 
a neutrino oscillation disappearance experiment using a single neutrino flavor in vacuum 
whereas in previous works p, 0] the method required the use of two different disappearance 
channels close to the first oscillation maximum. The possibility that the neutrino survival 
probability in vacuum can be sensitive to the neutrino mass hierarchy was first suggested 
in using reactor neutrinos in the context of high-Am 2 large mixing angle (LMA) solar 
neutrino solution. It was also pursued recently in the Hanohano project using a Fourier 
transform technique [l2| • We present here a complete formulation and physics understanding 
of the new method to determine the neutrino mass hierarchy. 

This method utilizes the feature that the atmospheric-Am 2 scale oscillation "interferes" 
with the solar-Am 2 scale oscillation in a different way depending upon the mass hierarchy. 
We first carefully re-formulate the difference as an advancement or retardation of the phase 
of the atmospheric-scale neutrino oscillation relative to the solar-scale oscillation; For the 
u e and i> e disappearance probabilities, the atmospheric oscillation is continuously advanced 
(retarded) for the normal (inverted) hierarchy. Since the core of the problem exists in 
extracting long wavelength solar oscillation out of the two short wavelength atmospheric 
scale oscillations due to Am| 2 and Am^, it is the critical matter to identify which quantity 
to be held fixed in defining the mass hierarchy reversal from the normal to the inverted, or 
vice versa. To our understanding, the only known appropriate quantity to hold is the effective 
neutrino mass squared difference Am 2 e (see Eq.fll])), the one which is actually observed in v e 
(or u e ) disappearance measurement at the first few oscillation maxima and/or minima 
We find it important to have the proper formalism because otherwise the sensitivity to the 
mass hierarchy resolution can be significantly overestimated. 

Then, we point out that an experiment using the resonant v~ absorption reaction 13| 
enhanced by Mossbauer effect recently proposed by Raghavan |l4[ gives the first concrete 
realization of the idea. (See [l5| for earlier proposals of Mossbauer enhanced experiments.) 
Because of the monochromaticity of the v e beam from the bound state beta decay (l6| . 
the phase advancement or retardation of atmospheric-scale neutrino oscillation is detectable 
after 20 or more oscillations if the source and the target are made sufficiently compact in size. 
The unique feature of this method is the use of this ultra-monochromatic neutrino source 
thereby sidestepping the very high energy resolution requirement requested with reactor or 
accelerator neutrinos. 

A concrete plan for the experiment may be as follows; In the first stage the atmospheric 
Am 2 e must be accurately measured in a way described in TtJ by a measurement at ~10 m 
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baseline. As we will see, the precision of the measurement must reach sub-percent level in 
order for this method to work. Then, in the second stage relative phase advancement or 
retardation of the atmospheric-scale neutrino oscillation must be detected by moving the 
detector to ~350 m, this is just beyond the solar-scale oscillation distance for 18.6 keV 
neutrinos. Based on this plan we give a detailed estimate of the sensitivity to the resolution 
of the mass hierarchy. 

Upon finishing the above measurement, it is natural to think about the continuation of 
the experiment by moving the detector to distances tuned for precision measurement of the 
oscillation parameters relevant for solar neutrinos, Am^ and #i 2 . We also give a detailed 
estimate of the accuracy for the determination of these parameters. It is notable that a 
resonant v e absorption experiment could allow to make precise measurements of all mixing 
parameters that appear in the v e — > u e oscillation channel: 6* 13 , # 12 , Am^, the effective 
atmospheric Am^, and the neutrino mass hierarchy. 

In Sec. [Ill we first describe the essence of the method in an explanatory way. We then 
define our theoretical machinery using the effective Am^,. We also clarify the basic require- 
ment for the measurement. In Sec. IIII[ we discuss the concrete set up of the experiment 
and describe our statistical method for the analysis. We then give the analysis results on 
sensitivity to the mass hierarchy resolution. In Sec. IIV} we describe the analysis for the de- 
termination of the solar oscillation parameters. In Sec. |Vj we finish by giving our concluding 
remarks. In Appendix |Al we give a different interpretation of the phenomena in terms of 
the effective mass squared differences whereas in Appendix |Bl we discuss how this method 
can be applied to the disappearance channel. This possibility has been discussed in 
however the baseline required is L ~ 11000 (-E/0.6 GeV) km. 



II. HOW CAN u e DISAPPEARANCE MEASUREMENT ALONE SIGNAL THE 
NEUTRINO MASS HIERARCHY? 

The vacuum v e survival probability using Ay = AmfjL/AE (Am-j = mf — m|) as 
shorthand for the kinematical phase, can be written without any approximation as 

P{v e ^v e ) = l-4|f/ e3 | 2 |t/ e i| 2 sin 2 A3i-4|f/ e3 | 2 |[/ e2 | 2 sin 2 A3 2 -P , 

= l-^sin 2 2# 13 [1- (c 2 2 cos2A 31 + S 2 2 cos2A 32 )] - P , (1) 

where the usual parameterization of the mixing matrix [19( is used with the abbreviations 
Cij = cos 6 \j and s y - = sin 6^-. Here E is the neutrino energy and L is the baseline distance. 
P Q , as defined by, 

P Q = 4|f/ e2 | 2 |[/ el | 2 sin 2 A 21 = sin 2 2# 12 4 sin 2 A 21 , (2) 

is the oscillation probability associated with the solar Am 2 scale. In this paper, we assume 
the standard three neutrino framework with CPT conservation (except for where we discuss 
possible CPT violation), which implies that Eq. (TTJ is also valid for v e . 

The neutrino mass hierarchy can be determined if one knows which kinematical phase, 
|A 31 | or |A 32 |, is larger. Namely, if |A 31 | > |A 32 | the normal hierarchy is the case, and if 
|A 3 i| < |A 32 | the inverted one. How we can know which is larger? When two waves with 
slight different frequencies travel together there arises the phenomenon of beating. The 
amplitude of superposed waves beats with the frequency corresponding to the difference of 
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two frequencies. For v e disappearance probability, Eq. flTJ), we have beating because the 
frequency of the two waves in ([1]) with |A 3 i| and |A 32 | differ by the frequency associated 
with the IA21I characterized by the solar Am^ scale, which is smaller by a factor of ~ 30 
compared to the atmospheric |A 31 | scale. One can decompose the superposed waves into the 
beating low frequency wave and the high frequency wiggles within the beat. Here, unlike 
the case addressed in most textbooks of basic physics, the amplitude of the two frequencies 
is different. Therefore, to tell which of the high frequencies has a larger amplitude, one 
has to determine whether the phase of the high frequency wiggles advances or retards as 
the wave evolves. This phase advancement or retardation can tell us whether the neutrino 
mass hierarchy is normal or inverted using only the disappearance measurement in a single 
channel, v e — > v e in vacuum. 

To demonstrate that it is possible, we plot in the left panel of Fig. [1] the oscillation 
probability P[y e — > z/ e ) as a function of the distance L from a source. As one can clearly see 
in the figure, the case of normal hierarchy (Am^ > 0, blue curve) starts to be distinguished 
from the case of inverted hierarchy (Am^ < 0, red curve) at around L ~ 300 m because 
of the advancement of the normal hierarchy's high frequency wiggles and the retardation of 
the inverted hierarchy's high frequency wiggles, i.e. the phase of the high frequency wiggles. 
What we mean by "beating" may be clearer in the left panel of Fig. [3] where the solar term 
is switched off. 

To see the difference between the oscillation probabilities of the normal and the inverted 
hierarchies more clearly, we define the quantity 



AP(u P 



PU 



u e ; Am 2 31 > 0) - P[v e 



v e \ Amjj < 0). 



(3) 



Where the hierarchy flip is defined so that |AP| is minimized for short baseline measure- 
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FIG. 1: In the left panel, antineutrino survival probability P{p e — ► Ve) is plotted as a function of 
L. The blue and the red curves are for the normal and the inverted mass hierarchies, respectively. 
In the right panel, AP/ \J P{B\z = 0), which is defined in the text, is shown. Here the neutrino 
energy was fixed to 18.6 keV. This choice comes from the experiment we will discuss later in this 
paper. 
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ments. (See Sec. Ill Al for the precise definition of the hierarchy flip.) We normalize AP by 
dividing by \J P(6i3 = 0) since -P(6>i 3 = 0) is the "background" for seeing #13 effects in v e 
disappearance. The quantity AP/ /P(6i3 = 0) is plotted in the right panel of Fig. [TJ The 
envelope of \AP\/y/P(6 13 = 0) slowly rises from to a plateau value of sin 2 26*13 at a distance 
about 25% further from the source than the first solar oscillation maximum (A21 = tt/2 oc- 
curs when L = 280 m) and remains reasonably constant at this value beyond the first solar 
oscillation minimum (A 2 i = tt). Clearly, even though the difference between the hierarchies 
is very small at short distances it becomes as large as theoretically possible, sin 2 2# 13 (see 
below), at large distances. 



Theoretical framework 



To decompose the superposed waves into the beat and wiggle, one has to correctly define 
the high frequency mode. We choose to do this by defining a short distance effective Am 2 
in the v e disappearance channel, see It is given by 



Am 



cos #12 1 Am 



'31I 



sin # 12 |Am 



'321 



(4) 



This is the only linear combination of Am^ and Am| 2 that can be determined by a short 
distance (over the first few atmospheric oscillations) v e disappearance experiment. Am 2 e is 
defined to be positive for both hierarchies. The choice of the atmospheric Am 2 is crucial; 
Extraction of the low frequency mode can be done correctly only if the high frequency mode 
is identified in an appropriate way. See Sec. Ill CI for more about it. 

Now, we rewrite P{y e —>■ v e ) without any approximation, using the positive and measur- 
able variables A ee = Am\ e L/\.E and A 2 i and the phase <fi & (which will be given below), as 
follows: 



P(v P 



sin 2 29 



13 



1 - \J 1 - sin 2 2#i2 sin 2 A 2i cos(2A ee ± 



P<s» (5) 



where the positive (negative) sign in front of the phase <fi & corresponds to the normal (in- 
verted) mass hierarchy; If Am 2 e is precisely determined, the only effect of the hierarchy 
manifest itself as this sign. Clearly, the amplitude of the sin 2 26 \% oscillations is modulated 
by y/l - sin 2 29 12 sin 2 A 2 i which varies from 1 to cos26>i 2 as the solar kinematic phase, A 2 i, 
changes from (mod7r) to tt/2 (mod7r). This amplitude modulation is the beat discussed in 
the previous section. 

The phase <p & in Eq. (J5]) depends only on solar parameters, A 2 i and #12, and given as, 



sin 1 



COS <p 



c 2 2 sin(2s 2 2 A 2 i) - s\ 2 sin(2c 2 2 A 2 i) 

y/l - sin 2 26»i2 sin 2 A 2 i 
c 2 2 cos(2s 2 2 A 2 i) + s\ 2 cos(2c 2 2 A 2 i) 
v 1 — sin 26 12 sin 2 A 2 i 



Or, alternatively 



arctan(cos 2Qyi tan A 2 i) — A 2 i cos 26*12. 



(6) 



(7) 



However, Eq. ([6]) contains quadrant information not available in Eq. (J7|). In Fig. (2], we show 
4> & as a function of A 2 i. Note that it is an odd, monotonically increasing function of A 2 i. 
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FIG. 2: (a) 4>& as function of A21. 0© = vr/2 when A21 ~ 1.9. (b) as function of L for various 
choices of the solar parameters within the current allowed region, ~ tt/2 at about 350 m 
somewhat independent of the precise values for the solar parameters. 



Due to the monotonically increasing nature of Q and the ± sign in front of , in Eq. (JSj) , 
for the normal hierarchy the high frequency wiggles undergo a phase advancement as the 
wave evolves, whereas for the inverted hierarchy there is a phase retardation 1 . It is easy to 
show that 



3 (A 2 i + 7r) =0 (A 21 ) + 27rsin 2 # 12 , (8) 



'0 



i.e. the phase advancement (normal) or retardation (inverted) is 2tt sin 2 #i 2 for every 7r 
increase of A 2 i- Eqs. (jSJ) and are the foundations of our investigation. 
Some important observations are worth emphasizing immediately: 



• OnlyA ee varies at the atmospheric scale. Everything else, including phase , varies 
at the solar scale. This is a useful distinction because these two scales differ by a factor 
of ~ 30. 

• The difference between probabilities corresponding to the two hierarchies ([3]) is given 
by 

AP 

, = = sin 2 2^13 sin 2 A ee sin , (9) 

V "(013 = 0) 

to leading order in #13. AP becomes visible when the phase develops to order unity. 
From Fig. [2] this occurs at around the first solar oscillation maximum, (A 2 i = 7r/2), in 
agreement with the features exhibited in Fig. [TJ From we also see that the normal 



Equivalcntly, one could interpret this phenomena as a change in the instantaneous effective Am^ tm . This 
interpretation is explored in Appendix A. 
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and the inverted hierarchies are maximally out of phase when 0© = n/2. This occurs 
when 

cot A 2 i cot(A 2 i cos 26»i 2 ) = - cos 20 12 . (10) 

The approximate solution to this equation, for 612 around 30 degrees, is A 2 i ~ § + 
cos 26' 12 / (vr sin 2 12 ). For sin 2 6*12 = 0.31, Eq. (TTOT) is satisfied when A 21 « 1.9, i.e., just 
beyond the first solar oscillation maximum. 

• We notice that by doing transformation 6 12 — > n/2 — 9 12 (or s 12 <-> c 12 ) in flB]) flips 
the sign of </> . It means that the mass hierarchy is confused by the presence of the 
particular type of CPT violation, the subtle one discussed in (2o| . It occurs because we 
use the antineutrino channel; If we are able to use the neutrino channel, the problem 
of confusion does not arise because we know that B\2 < 7r /4 from the solar neutrino 
observation. 



B. Experimental strategy and requirements to the accuracy of measurement 

We now describe our experimental strategy for determining the hierarchy. First one needs 
to make a very precise measurement of Am 2 e at short distances from the source, at around 
half the atmospheric oscillation length, Lq/2 where 



L = 



4ttE 
Am 2 ce 



18.4 



E 



18.6 keV 



Am 2 



2.5 x 10" 3 eV 2 



m . 



(11) 



Then, one must determine the sign in front of where the phase difference between normal 
and inverted hierarchy is maximal, i.e. A 2 i ~ 2 as can be seen in Fig. [2j With A ee /A 2 i ~ 30, 
it can be translated into the distance L ~ [(2 x 30)/tt]Lq ~ 20 Lq away from the source. 

How well do we need to know Am 2 e ? If we do not know the value of Am 2 c with sufficient 
precision, one can adjust Am 2 e , depending upon the normal (NH) and the inverted (IH) 
hierarchies, so that the relation 



2A, 



ee NH 



+ 



2A 



cc IH 



(12) 



holds at some L/E. Then, from Eq. (jSJ) the hierarchy is completely confused as we will 
explicitly see in Sec. IIII Al It can be translated into 



Am 2 



1 111 



— Ami 



NH 



Am 2 , 



'0 



(13) 



which leads to the requirement for the accuracy of measurement of Am 2 c as 
<5(Am 2 e )/Am 2 e < 2Lq/-kL ~ 3%. This is the critical value at which the hierarchy is com- 
pletely confused, and if we want to do better, e.g., at least a factor of 3 times better than 
3%, the first requirement is ~1% measurement of Am 2 e . 

At this point the astute reader will be concerned that we are attempting to determine a 
oscillation phase after 20 or so oscillations. Aren't the atmospheric oscillations completely 
washed out long before 20 L ? To observe atmospheric oscillations 20 or so oscillations 
out, assuming coherence, demands that we can determine L/E per event to high precision. 
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FIG. 3: The averaged P{v e — ► f e ) versus distance measured in number of oscillation lengths, Lq. 
To help with visualization Pq = sin 2 29u cos 4 #13 sin 2 A21 has been set to zero for these figures, (a) 
0.1% Gaussian resolution on L/E. (b) 0.5% resolution on L/E. (c) 2% resolution on L/E. In (b) 
the oscillations are reduced at 20 or so oscillations out, but still observable. Whereas in (c) the 
oscillations are averaged out by 15 Lq. 



Figure [3] shows how the oscillations are reduced as the resolution on L/E is increased and 
that a precision better than 1% is required to see oscillation at ~ 20 oscillation lengths from 
the source. 

In summary, to measure the sign of the phase and determine the hierarchy we have two 
important initial requirements: 

• The resolution on Am 2 e measured in a short baseline experiment (~ L /2 from the 
source) must be ~1% or less. 

• The resolution on L/E for the long baseline experiment (~ 20 L Q from the source) 
must be ~1% or less. 

Usually, the neutrino energy cannot be determined experimentally with this accuracy, there- 
fore, these preconditions are very difficult to meet, if not impossible. 



C. Mass hierarchy reversal and comments on the reactor neutrino method 

We have emphasized in the foregoing discussions the importance of identifying the quan- 
tity to be held fixed to define the mass hierarchy reversal, and proposed Am 2 e as the solution. 
In this subsection we want to elaborate this point and clarify how the difference between 
the normal and the inverted hierarchies can be made artificially larger by choosing different 
variables to define the hierarchy reversal. 

In Fig. HI presented is the survival probability of electron anti-neutrinos at a baseline 
of 50 km from a source which is averaged over the uncertainty of energy 3%/v/E/MeV as 
a function of neutrino energy E. It is shown in the left panel of Fig. H] that if we hold 
|Am| 2 | fixed in reversing the hierarchy from the normal to the inverted (as was done in [lll |) 
the difference between the normal and the inverted hierarchies is clearly visible. However, 
the obvious distinction seen in the left panel of Fig. H] disappears if we use Am 2 e to hold 
when switching between the hierarchies as shown in the middle panel of Fig. HI A simple 
explanation for such a marked difference is that by holding |Am 2 2 | fixed in reversing the 
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hierarchy from the normal to the inverted one maps the neutrino mass spectrum into a 
significantly different one, lAm^-Jin = jAm^lNH ~~ 2Am| 1 . (See Appendix lAl for more about 
it.) This artificially enhances the difference between the hierarchies, as demonstrated in 
Fig. H 
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FIG. 4: The z/ e disappearance probability at 50 km from a reactor source as a function of the 
neutrino energy. The resolution on the energy is assumed to be 3% / y / E /MeV . The blue (light gray) 
and the red (dark gray) curves are for the normal and the inverted mass hierarchies, respectively. 
In the left panel |Am| 2 | is held fixed when the hierarchy is flipped, whereas in the middle panel 
A?7ig e is held fixed. In the right panel Ara^ for the inverted hierarchy is 0.6% larger then for the 
normal hierarchy. 



Despite that there exist some discrepancies at energies around 4 MeV, a peak energy of 
signal events in reactor experiments, the two disappearance probabilities can be made very 
close to each other over the energy range of 4 to 6 MeV if we allow the value of Am^, for 
the inverted hierarchy to increase by 0.6% relative to the value for the normal hierarchy, as 
can be seen in the right panel of Fig. |H The lesson we have learned from this exercise is 
that unless we know Am^ e to a sub-percent level, the hierarchy can be easily confused by a 
tiny shift of the atmospheric Am 2 . 

One might want to use neutrinos from reactors to determine the mass hierarchy by the 



present method, as suggested by the authors of 11] . However, the experimental require 



ments are very severe in order for this method to work. We believe that the three points 
discussed above deserve attention if one want to make a reliable estimate of sensitivity to 
the mass hierarchy resolution with reactor neutrinos; (1) definition of hierarchy reversal, (2) 
implementation of realistic energy resolution, and (3) robustness of the sensitivity to mass 
hierarchy resolution against the uncertainty of Am^. 

The recent Hanohano proposal, [l2|, sidesteps some of these requirements by using a 
sophisticated Fourier transform technique. Detailed analysis around the atmospheric peak in 
Fourier space can tells us whether |Am3 2 | is greater or smaller thanlAm^l, thus determining 
the hierarchy. There is no need to explicitly deal with the definition of hierarchy reversal. 
In this way, the two curves shown in the two right panels of Fig. H] can, in principle, be 
distinguished. If this can be achieved in a realistic setting of the experiment, it represents 
a remarkable advance for reactor neutrino experiments, demonstrating the power of the 
Fourier transform technique. 
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III. EXPERIMENTAL SETUP, ANALYSIS METHOD AND RESULTS 



Most probably, the unique possibility to satisfy both requirements of the previous sec- 
tion is to make use of monochromatic beam from the bound-state beta decay 3 H — > 
3 He + orbital e~ + v e and detect it by the time-reversed resonant absorption reaction 
v e + 3 He + orbital e~ — ► 3 H, the idea recently revived by Raghavan [3]. For the Mossbauer 
enhancement to work the method requires AE/E better than 10 -16 , a truly monochromatic 
neutrino beam of 18.6 keV. Thus, the energy spread of the beam in such an experiment is 
automatically much narrower than what is required. 

Having monochromaticity satisfied, the remaining requirement is on the baseline L. The 
short and long baseline measurement are to be done at the distances of ~10 m and 20Lo ~ 
350 m from the source, respectively. Then, the requirement on AL/L of, say, 0.3% implies 
that the size of the source and the target must be limited respectively to ~3 cm and ~1 
m in short and long baseline measurement. The former places severe limits on the sizes of 
the source and the detector. Assuming this requirement can be met, it has been shown by 
analyzes with several concrete settings that a sensitivity to Am^ of ~ 0.3 (sin 2 29i 3 /0.1)~ 1 % 
at la CL is possible [13]. Therefore, the first requirement listed in Sec. IHBI can be met if 
sin 2 2^13 > 0.05. 



A. Experimental setup and how two-stage measurement works 

We consider a possible experimental set up for determination of the neutrino mass hier- 
archy based on a Mossbauer enhanced u e absorption experiment [14|. By using the revised 
estimate of the cross section a res = 0.3 x 10~ 32 cm 2 [ijj the event rate is given by 



(S M \ ( L \ ~ 2 
-^vtt- — - — day" 1 . (14) 

lMCi- 100 gj \10mJ J y J 

Therefore, one obtains about one million (2500) events per 3 days (10 days) for 1 MCi source 
and 100 g 3 He target at a baseline distance L = 10 (350) m. The experiment should be done 
in two stages; In the first stage, Am 2 e and sin 2 2#i3 should be determined very precisely by 
a measurement at baseline around ~10 m as demonstrated in 17 . We thus do not discuss 



it further by just assuming the setting with 10 location measurement called the run IIB in 



17| . For the second stage, we examine the setting with detectors placed at five locations 
L = 350 ± (0, 5, 10) m from the source, which covers a complete atmospheric oscillation 
length. The results of these two stages should be combined as we will discuss bellow. 

In Fig. we present the expected number of events for Am 2 e = 2.5 x 10 -3 eV 2 and 
sin 2 2^13 = 0.1 at 5 locations mentioned above for the normal hierarchy, for 2000 events 
at each detector location in the absence of oscillation, which correspond to an exposure 
of about 20 MCi-g-day. 2 Throughout our analysis except for that in Sec. HVl we fix the 
solar parameters to their best fit values as sin 2 #i 2 = 0.31 and Am^ = 8.0 x 10~ 5 eV 2 



2 In view of Fig. [5) one may think that measurement at least 3 different baseline distances would be enough 
to detect relative phases of the beat and the wiggle around it. Though it is true, it turned out that the 
five point measurement can achieve better accuracy than the three point measurement for an equal total 
number of events. 
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and do not try to vary them. Am^ e is taken with a tentative expectations for its accuracy, 
Am 2 c = 2.5 x 10~ 3 (1 ± 0.5 %) eV 2 . The case of the normal and the inverted hierarchy are 
plotted by the blue (light gray) and the red (dark gray) curves, respectively. It appears 
that, if enough statistics is accumulated and Am 2 e is separately measured accurately, the 
five point measurement, can discriminate between the normal and the inverted hierarchies. 
The details of the statistical procedure is described in the next subsection. 

We show in three panels in Fig. [6] the allowed regions in the sin 2 2#i 3 — Am 2 e space 
obtained by the first stage alone (upper panels), by the second stage alone (middle panels), 
and the one obtained when the first and the second stages are combined (lower panels). The 
input parameters are taken as sin 2 2# 13 = 0.1, Am 2 e = 2.5 x 10~ 3 eV 2 , and the mass hierarchy 
is assumed to be the normal one. In the first-stage measurement to be done at around the 
distance Lq/2 ~ 10 m from the source, sin 2 2# 13 as well as Am 2 c can be determined rather 
precisely. But, there is no sensitivity to the mass hierarchy because the baseline is too short, 
as is the case of reactor # 13 experiment. 

The feature of the second-stage experiment presented in the middle panels of Fig. [H] 
requires some explanations. First of all, there exist many isolated island-like allowed regions. 

sin 2 0, 2 = 0.31 Am", = 8.0 x 10~ 5 eV 2 sin 2 20, 3 = 0.1 
650 
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FIG. 5: Plotted are the expected number of events to be collected by detectors placed at the 
distances 340, 345, 350, 355 and 360 m from the source for Am 2 e = 2.5 x 10~ 3 eV 2 and sin 2 26>i 3 = 
0.1 for the normal hierarchy, indicated by solid circles with error bars. Here Pq is included (solar 
parameters are fixed to their best fit values). Expected number of events in the absence of oscillation 
is assumed to be 2000 at each detector position. The blue (light gray) solid curve which passes 
through the data points indicates a theoretical expectation assuming the normal hierarchy with 
Am 2 e = 2.5 x 10 -3 eV 2 and sin 2 2#i 3 = 0.1 whereas the red (dark gray) curve off the data points 
is the one of inverted hierarchy with same values of Am 2 e and #i 3 . Assuming that Am 2 e is known 
to ±0.5% precision, the thick-dotted (thin-dotted) lines above and below the solid lines are the 
expectations with Am 2 e 0.5% above (below) that value. 
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FIG. 6: The allowed regions in the sin 2 26±3 — Am e 2 e plane obtained by the first stage alone (upper 
panels), the second stage alone (middle panels) and the combined result of these two stages (lower 
panels) for the case where the input parameters are sin 2 2#i3 = 0.1, Am 2 e = 2.5 x 10~ 3 eV 2 with 
the normal hierarchy, based on the x 2 analysis described in Sec. IIIIB1 The solar parameters are 
fixed to their current best fit values. The first stage is assumed to be performed at 10 positions at 
around L ~ 9 m following the run IIB in 17J whereas the second stage is assumed to be performed 
at 5 locations at L = 350 ± (0, 5, 10) m from the source. The three closed curves (red, blue, green) 
from inner to outer denote the ones obtained at 90% CL, 95% CL, and 99 % CL. 



Suppose that #13 is given from some other experiments. Then, due to the cosine term in 
Eq. ((Sj) the measurement will give a periodic solution 

Am^Am c 2 c + ^, (n = ±1, ±2, ...) (15) 

for both of the hierarchies. If we take approximation Q = tt/2, there is a shift 2ttE/L 
between the adjacent Am 2 e solutions for the normal and the inverted hierarchies. Hence, the 
solutions of Am 2 e are alternating in the normal and the inverted hierarchies with a constant 
shift 2itE/L between the adjacent solutions, the feature clearly visible in the middle panels 
of Fig. The contours have finite widths, though prolonged, along the sin 2 2t9 13 direction 
because the second stage measurement alone can determine t9 13 but only with a limited 
accuracy. 

It's clear that the second stage alone can not determine neither the value of Am 2 e nor the 
mass hierarchy. However, once the results of the first and the second stages are combined 
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all the fake solutions are eliminated, as demonstrated in the lower panels of Fig. [61 It is 
evident that Am^ e must be determined with high accuracy in the first stage measurement. 



B. Analysis method and definition of y 2 



Now we give details of our analysis method. To calculate the sensitivity of the deter- 
mination of the neutrino mass hierarchy by the first and second stages of the experiment 
combined, we compute 

A Xmin( sin2 2 #i3) = Xmin( wron g hierarchy) - Xmin( true hierarchy) , (16) 
where Xmm(t rue / wron § hierarchy) is the minimum of 

X 2 (true/ wrong hierarchy) = \\ + xl(^ Iue I 'wrong hierarchy) , (17) 
with the first stage \ 2 given by 



Xi s t(AmL sin 2 29 
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(18) 



and the second stage \ 2 given by 



X2nd 



(Am 2 ee , sin 2 2^ 13 ,sign(Am 2 1 )) = £ 
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(19) 

' "° s is the number of events at baseline Lj for a given input hierarchy and values of the 
mixing parameters (observed events), while N^ ,7rt eo is the theoretically expected number of 
events at baseline Lj for a given set of oscillation parameters. We use, as in Ref. 
correlation matrices 



s 



a 



1 



(7, 



[i + (E* 4>. 



3j, the 
(20) 



where a 2 , = a 2 sys + 1/AT 



1,2-th 



The estimation of the uncorrelated systematic uncertainty is the key issue in our sensi- 
tivity analysis. In [17] two cases of the error, cr usys = 0.2 % and 1%, are examined. The 
former corresponds to the case of movable detector, assuming that a direct extraction or 
real-time detection of produced 3 H atom is possible. The latter is a tentative estimate for the 
case of multiple identical detectors. We use the correlated systematic uncertainty a c = 10 
%, though the results are insensitive to a c . We adopt the same numbers as these for the 
uncertainties in our analysis for both the first and the second stages of the experiment. 

In our x 2 analysis, for simplicity, the solar mixing parameters are fixed to the current 
best fit values, ignoring their uncertainties. Therefore, the fitting parameters are: Am 2 e , 



sin 29 13 and the sign of Am^. Note, however, that x 2 s t f° r the first stage has practically 



no dependence on solar parameters as well as the mass hierarchy due to too short baseline. 
The current uncertainty on the values of the solar parameter are not so relevant for the 
determination for the hierarchy because <j) & does not depend strongly on them (see Fig. [2b) 
as pointed out previously. 
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C. Sensitivity to mass hierarchy determination 
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FIG. 7: The region of sensitivity to resolving the mass hierarchy in sin 2# 13 — event number (per 
detector) space. The black solid, the red dashed, and the blue dotted curves denote the region 
boundary at 90%, 95%, and 99% CL, respectively. The left and right panels are for the cases of 
uncorrelated systematic uncertainties of 0.2% and 1%, respectively. 



Now, we present our estimation of the sensitivity of the two-stage Mossbauer v e experi- 
ment to the determination of the neutrino mass hierarchy. In Fig. [7] we present the region of 
sensitivity to resolving the mass hierarchy at 90%, 95%, and 99% CL, denoted by the black 
solid, the red dashed, and the blue dotted curves, respectively. The abscissa is sin 2 2# 13 and 
the ordinate is the number of events per detector in long baseline experiment. Here, 90%, 
95%, and 99% CL correspond to A^ 2 ^ = Xmi n ( wron g hierarchy) — ^^(true hierarchy) = 
2.71, 3.84 and 6.63, respectively. Left and right panels are for the cases where the uncorre- 
lated systematic uncertainty is 0.2% and 1%, respectively. 



Assuming a 



usys 



0.2 %, we can see that the hierarchy can be determined at 90 (99) 
CL for sin 2 29 13 down to 0.035 (0.06) for 1000 events and down to 0.025 (0.04) for 2000 
events. For a larger error a usys = 1%, the critical value of sin 2 2#i 3 will be ~ 0.08 (0.11) at 
90 (99) % CL approximately independent of the number of events larger than ~ 1000. This 
is caused by the fact that the precision on the determination of Am 2 e at the first stage is 
highly dependent on this systematic uncertainty [17 . 



IV. DETERMINATION OF THE SOLAR PARAMETERS 

In the second stage of the Mossbauer experiment, the detectors are placed just after the 
first solar oscillation maximum in order to determine the mass hierarchy, as we discussed 
in the previous section. In a possible third stage of the experiment, one could envisage 
moving the detectors somewhat closer to the source in order to cover the region around this 
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maximum. It will allow us a precise determination of the solar-scale oscillation parameters, 
Anijj and On- 

In order to optimize the determination of these parameters, we assume that the measure- 
ments will be performed at the following 10 different detector locations, 

L n = [200 + 50 (n - 1)] m, (21a) 
L' n = L n -L /2, n = 1,2,. .,5. (21b) 

First, in order to observe the oscillation driven by the solar Am^, we consider the configu- 
ration of 5 detector locations (L n ) separated by 50 m, ranging from 200 m to 400 m as in Eq. 
(1 2 laft . In this way, we can cover the whole range of solar-scale oscillation before and after the 
dip due to the oscillation maximum (see Fig. [1]). Second, in order to minimize the unwanted 
effect due to the rapid oscillations driven by the atmospheric Am 2 e , we have to place, for 
each location L n , another detector (or move detector to another location) at L' n separated 
from L n by half the atmospheric oscillation length, L /2. See Eq. (121b|) . The setting of five 
pairs of detector locations works because the average of the oscillation probabilities corre- 
sponding to these two positions L n and L' n is, at first approximation, independent of 0\z an d 
Am 2 e and thus is free from any ambiguity associated with the atmospheric oscillations. 

In this calculation we assume that the exposure is such that the same number of events 
(1000) in the absence of oscillation are collected at each detector location. We take the 
pessimistic value of the systematic uncertainty cr usys = 1% instead of 0.2 % because the 
results do not differ much; The sensitivity is still dominated by the statistics in this regime. 
As in the previous sections, the input values of the other oscillation parameters are fixed as 
Am 2 ee = 2.5 x 10~ 3 eV 2 , sin 2 12 = 0.31 and Am 2 21 = 8.0 x 10~ 5 eV 2 . 

In Fig. [8] we show how well Am^j and sin 2 On can be determined by this experiment. 
The obtained sensitivity contours with sin 2 # 13 = are presented by the green lines (which 
practically overlapped with the black ones) in Fig. [HJ From this result, we conclude that the 
precision one can achieve for sin 2 On is about 1.6 (4.0) % at 68.27 (99.73) % CL for 2 DOF 
whereas that for Am^ is about 0.9 (2.2) % at 68.27 (99.73) % CL for 2 DOF. Comparing 
this result with the sensitivity which can be achieved by the dedicated reactor experiment 



for the similar number of events studied in Ref. [211 ]. the precision we obtained for sin 2 0i 



7 12 

and Am^ in this work are comparable or better, in particular, for Am 2 x (see e.g. Fig. 5 of 



the first article in 211 ]) 



In the above argument, we have assumed that L (or Am 2 c ) can be determined very 
precisely in the first stage of the experiment allowing us to place the detector pair separated 
exactly by Lq/2. In reality, however, the determination of Lq will always carry some uncer- 
tainties. Moreover, if sin 2 20i3 is too small or the uncorrelated systematic uncertainty (cr usys ) 
is not small enough, the neutrino mass hierarchy will not be determined in the second stage. 
Then, the question is; What is the additional ambiguities that occur in determination of the 
solar oscillation parameters in the presence of these uncertainties? Fortunately, we can show 
that our choice of detector locations defined in Eqs. (121 al) and (121b|) guarantees that these 



uncertainties will not deteriorate the precision in the determination of the solar parameters. 

In order to check the stability of the results against the uncertainties mentioned above, 
we performed a fit by assuming the incorrect mass hierarchy and/or slightly incorrect value 
of Am 2 c by shifting its value by 4%. The results are shown by the blue, the magenta and 
the light blue curves in Fig. [8] where we can see that for each CL, all the curves are very 
close to each other, indicating that our results are quite stable against these uncertainties. 
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FIG. 8: Expected allowed regions in the plane sin 2 612 — Am 2 ^ for 1000 events taken at each 
detector position at 68.27% (inner curves), 95.0% (middle curves) and 99.73% (outer curves) CL 
for 2 DOF. We assumed the measurement will be done at 10 different detector locations, L n = 
[200 + 50 (n — 1)] m and L' n = L n — Lq/2 with n = 1,2, ...,5. The input values used are: 
sin 2 #12 = 0.31, Am^ = 8.0 x 10~ 5 eV 2 , Am 2 c = 2.5 x 10~ 3 eV 2 . The green and black curves 
are for the input values of sin 2 2#i3 = and 0.05, respectively. The other curves are also for 
sin 2 2#i3 = 0.05 but fitted with: the wrong hierarchy (red), Am 2 e shifted by 4% (magenta) and the 
wrong hierarchy and A?n 2 e shifted by 4% (cyan). Here we have supposed crusys — 1%. 
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We also examine the stability of the sensitivities to possible nonzero values of #13. In 
Fig. Owe show the sensitivity contours for the input values of sin 2 2# 13 = 0.05 (black curve) 
in comparison to that with sin 2 2#i3 = 0.0 (green curve). We note that these 2 curves 
are very close to each other implying that the effect of nonzero #13 is quite small. We 
emphasize that this feature is in a marked contrast with any of the alternative methods 
for measuring #12, including the reactor 21] and the solar neutrino methods 22|. In these 



methods there is no way, unless determined by some other experiments, to control the 



errors which comes from nonzero #13 [23j, whereas in our present method, there is a built-in 
mechanism which eliminates the uncertainties that come from the atmospheric oscillation 
parameters, in particular from #13. 

In Fig. Owe show how the precision for 1 DOF of each of the solar parameters improves 
with statistics, that is, as a function of the number of events that would be collected at each 
detector in the absence of oscillation. As in the previous plot, we have assumed that cr usys = 
1%. We observe that until ~ 1000 events, the precision improves rapidly as the number 
of events increase and then after, the rapidness of the improvement become smaller. The 
ultimate precision which can be achieved in the limit of very large statistics with ~ 10000 
events (per detector) for cr usys = 1% would be ~ 0.4 (1.2) % for 68.27% (99.73%) CL for 
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exposure at 68.27% and 99.73% CL for 1 DOF. Here <r usys = 1%. The dashed lines show the input 
values. 



sin 2 #12 and ~ 0.2 (0.7) % for 68.27% (99.73%) CL for Am 2 21 . Notice that measurement 
of Am 2 to the accuracy of sub % level is usually not possible because the absolute energy 
scale of leptons is hard to be determined at such accuracy, and it acts as a limiting factor 
for reducing the Am 2 error [7]. In our case, monochromaticity of the v e beam with known 
energy removes the necessity for such calibration for determining absolute energy scale. 



V. CONCLUDING REMARKS 

In this paper, we have shown that the neutrino mass hierarchy can be determined only 
by a single experiment based on resonant v e absorption reaction enhanced by the Mossbauer 
effect. In its first stage, Am 2 e must be accurately determined by L ~10 m measurements. In 
the second stage of the experiment, the relative phase of the atmospheric oscillation (wiggles) 
between the normal vs. inverted mass hierarchies, phase advancement or retardation must 
be discriminated by measurements at L ~350 m. It can be executed thanks to the ultra 
monochromatic nature of neutrino beam which would allow us to detect relative phase of 
atmospheric oscillation even after 20 oscillations. We stress that the monochromaticity is 
special to the resonant v e absorption reaction, most probably the unique experimental way 
to realize this principle. 

We concluded that assuming cr usys = 0.2 %, the hierarchy can be determined at 90 (99) 
% CL for sm 2 2# 13 down to 0.035 (0.06) for 1000 events collected at five different positions 
and down to 0.025 (0.04) for 2000 events. For a larger error a usys = 1%, the critical value of 
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sin 2 2^13 will be ~ 0.08 (0.11) at 90 (99) % CL approximately independent of the number of 
events for more than ~ 1000 events. 

We have also shown that in a third stage of the experiment, by positioning detectors in 
different locations, it is possible to determine the solar-scale parameters too. The precision 
one can achieve with 1000 events in each detector for sin 2 # 2 i is about 1.6 (4.0) % at 68.27 
(99.73) % CL whereas that for Am 2 21 is about 0.9 (2.2) % at 68.27 (99.73) % CL, all for 
2 DOF. The accuracies improve as data taking proceeds and the errors would reduce by a 
factor of ~ 3 when an order of magnitude larger number of events is taken. 

It is in fact remarkable that the experimental set up considered in this paper allows us to 
make a precise measurement of all the mixing parameters appear in the oscillation probability 
for v e —>■ v e channel, #12, #13, Am^, Am 2 e and the neutrino mass hierarchy only by a single 
experiment, independent of any other experiments. The only parameters inaccessible by 
the method are #23 and the CP phase 5, but precision measurement of the other mixing 
parameters will certainly help in determining these parameters by other experiments. 

Some remarks are in order: 

• As noted in Sec. Ill A[ our method of determining the mass hierarchy can be confused 
by the CPT violation of the type with 612 in different octants in the neutrino and the 
antineutrino sectors. It in turn implies that once the mass hierarchy is determined 
by some other means, the measurement we propose in this paper can be regarded as 
testing that particular type of CPT violation. 

• Due to the capability of resolving atmospheric wiggles even after 20 atmospheric oscil- 
lations away from the source, this experiment could be sensitive to other dynamics 
which differ from the standard three flavor mass induced oscillation one. In particu- 
lar, one can perhaps probe mechanisms such as: neutrino decay [25| , decoherence [26| , 
mixing with sterile neutrinos [27j . violation of Lorentz invariance or the equivalence 



principle [28(, long-range leptonic forces [29J, extra-dimensions [3 



Finally, we emphasize that realizing the proposed Mossbauer v e experiment is crucial in 
order to make all the discussions in this paper meaningful. The principle of the experiment 
is well formulated by having the decay producing source v e and the absorption reaction in 
the target as T-conjugate channels with each other, allowing mutual cancellation of large 



part of the atomic shift 14], |3l| . Nonetheless, it still remains unclear to what extent the 3 He 
atom in the target experiences exactly the same environment as the 3 H atom in the source 
32| . This point has to be (and will be) investigated in a systematic way. If the experiment 
is indeed realizable it will certainly open a new window for neutrino experiments which can 
explore the neutrino properties in otherwise inaccessible manner. 
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APPENDIX A: THE EFFECTIVE ATMOSPHERIC Am 2 

The advancement or retardation of the phase, , could be interpreted as a change in 
the effective atmospheric Am 2 . Such an effective atmospheric Am 2 , at L/E, is naturally 
defined by 

a 2 {r /T?\ - d(2A cc ±(j) & ) 2 1 2 sin 2 26 12 sin 2 A 2 l v 

Am eS {L E) = rnoln = Am «= ± 5 Am 2i cos20 12 , 2 Al 

a{L/2E) 2 1 — sin 26 12 sin A 2 i 

4 - ~ sin 2 26»i2 sin 2 A 21 \ „ /s 2 2 - ~ sin 2 26» 12 sin 2 A 21 \ . . 2 . 

1 l Am 3il + — i r^~r |Am 32 |. 



1 - sin^ 2^12 shr A 21 J ' J1 1 \ 1 - sin 2 2^ i2 sin 2 A 21 j 

Where the positive (negative) sign, in Eq. (I All) , is for the normal (inverted) hierarchy. For 
A 2 i < 1 we recover Am 2 e , since this is in essence its definition, see Ref. [oj for details. Thus 
Eq. (IA.1I) is the generalization of Am 2 e to arbitrary L/E. In Fig. [10] we show the variation 
of Am 2 ff as a function of baseline for both the normal and inverted hierarchies. The extrema 
of the difference between Am 2 s (L/E) and Am 2 e are given by 

1 2 oa sin 2 2fli2 sin 2 A 21 JO at A 2 i = nn 

2 Am21 COS ^ 12 l-sin 2 2^sin 2 A 21 " \ \ gg? Am 2 , at A 21 = (2n - 1)tt/2. (A2) 

As the oscillations go through the solar maxima, the fractional increase (decrease) in the 
effective atmospheric Am 2 is 



1 sin 2 26»i2 Am^ 

2 cos 26*12 Am 2 



4% , (A3) 



for the normal (inverted) hierarchy. It is the accumulated effect of this small difference in 
the effective atmospheric Am 2 for the normal and inverted hierarchies that leads to the 
observable phase difference shown in Fig. [TJ Since 



2 

as it must. 



l - d(L/E) Am 2 ff (L/£) = (2A CC ± <f, Q )\ L/E , (A4) 



APPENDIX B: DISAPPEARANCE 

For U/j, disappearance in vacuum, most of the results of this paper can be applied, using 
the following translation 

sin 2 fl 12 =» jf^ 2 and sin 2 6 13 |^ 3 | 2 - (Bl) 
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FIG. 10: Am^g as function of baseline at 18.6 keV neutrino energy for both the normal and inverted 
hierarchies. The dashed line is Am^ which has been assumed to be the same for the normal and 
inverted hierarchies for this plot. In general, unless Am^ is determined precisely it need not take 
exactly the same value for both hierarchies. The maximum difference between Am^ and Am^ 
occurs at the solar oscillation maxima, i.e. when A21 = (2n — l)vr/2. 



With this translation the atmospheric Am^ and phase, , are given by 
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Since numerically 
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the phase for v„ disappearance is approximately minus the phase for v e disappearance, i.e. 



(B6) 
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Therefore the normal hierarchy is retarded while the inverted hierarchy is advanced. This 
is the opposite sign to what happens for v e disappearance but the rate of change is approx- 
imately the same. 

The effective atmospheric Am 2 for disappearance at any L/E is given by 



Am 2 ff (/z/z) 
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Again for small A21 this reduces to Am^ given by Eq. (1B2I) . 

For v T disappearance, just replace | | 2 with \U T A 2 in the equations of this appendix. 
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